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Earth Resources Survey Systeins are used to gather 
infornation to ticlp governrent, private or foreign 
users with statutory or other requirements to solve 
probicns or manage resources. These systems -consist 
of remote sensors in aircraft or spacecraft, data 
formatting and telemetry links, preprocessors, 
extractive processors, and user applications models. 
Remote sensors register attributes of the terrain 
being covered. The sensor data may be relayed to 
ground facilities through telemetry linUs. At the 
facilities, preprocessing and extractive processing 
arc performed, to extract information from data. 

Finally, user applications models are used to relate 
iaforciation outputs from the extractive processing, 
end ancillary data, to the user in tortus most iacaning“ 
fui to him. V 

Vith the- launch of the EKTS-1 satellite In July 
1972, the ability of renoto sensors to collect periodic 
tmltispcctral data was significantly increased over 
■ the previous and continuing airborne sensor capability. 
The potential of coverage of a given site on the 
earth every 18 days, facilitates change detection and 
opens up possibilities of analysis techniques using 
the time signature of remote sensible terrain attri- 
butes* Finally, the ability of the ERTS system to ^ 
synoptical ly view 100 n.ni^ frames'* (with 7.2 x lO'* 
temples per frame) facilitates surveys of large areas, 

= end because of this capability, more users of Eai;th 
•Resources System data arc actively interested. 

At the sane time however, considerable strain 
Is being put on the existing preprocessing and 
extractive processing hordvare and methodology because 
of the potential of ERTS for spcctral-spatial-temporal 
. analyses, because of the high data rates involved, 
because of the need to correlate data with existing 
airborne sensors for the high spatial resolution view 
often required, and because of user-irposed timeliness 
requirements. Additionally, considerable development 
of the user applications model area is required if 
the full potential of Earth Resources Systems to 
deliver cost-bencf icial results is to be achieved. 

At ERIM, we have been exploring for ten years the 
utility of renote sensing data to various uatuiui 
resources problems, to developing the system method- 
ology and hard'^ave necessary to rapidly provide 
accurate Information to users at low cost, and 
developing user applications i»jodels to enhance the 
value of Earth Resources Systcr^ data. Some recent 
developments in* extractive processing methodology and 
hardware and in user applications psodcl development 
arc the subject of this paper. 

Techniques for spcct ral-tcnporal-spat lal proces- 
sing were developed to permit rtorc reliable recogni- 
tion of daSGCS of vegetation and other terrain 
categories with distinctive tlmC-varying spectral- 
spnLlaJ signatures. Ejtamptcs of wetlands mapping and 
land resource inventory will be presented and discussed. 

To porfom extractive procecsing rapidly to meet 
user timeliness requ i rcBefi>ts , a high speed parallel 


digital special purpose processor MIDAS (Multivariate 
Interactive Digital Analysis System) is being developed 
at ERIM under R/\SA funding. This prototype processor 
will Implcncnt maximum likelihood class i f ica tion , with 
versatile data processing, at rates comparable to the 
data collection rate of ERTS. The design philosophy 
of this processor will be discussed in this paper. 

Last, an example of a user model developed to 
predict the yearly production, of mallard ducks , an 
important migratory waterfowl, from remote sensing 
and anclllory data, will be described. 

INTRODUCTION^ 

With the launch of the Earth Resources Technology 
Satellite (ERTS-1) in late July 1972, the attention 
of over 300 NASA-funded scientists focused on the 
application of data from the return beam vidicon (REV) 
and multispect ral scanner (Mf?S) to earth resources 
problems In the United • States and throughout the 
world. For those of us involved in the design of 
prototype earth resources systems to assist in resource 
management, ERTS provided a synoptic, hopefully 
periodic, broad scale look at the earth's resources 
attainable only with great difficulty from previously 
available aircraft sensors. 

While data from ^aircraft sensors will still be 
sought and used by resource managers because of Its 
relatively high spatial resolution, ERTS-likn data 
from satellites is definitely here to stay. ERTS, 
because ot its wide area coverage, has renewed the 
interest of resource managers in operational earth 
resource systems t^hich can supply them data on a 
routine basis. For systems designers , this is a 
challenge for three reasons: 1) The data -rate from 

earth resources sensor systems increased by over an 
order of magnitude with the launch of ERTS-1, creating 
problems of data acquisition and processing; 2) the 
demands for timely information from an earth resources 
system call for throughput requirements at least four 
orders of magnitude greater than can be obtained with 
currently implemented . algorithms on general-purpose 
digital computers 3) what is easily extracted from, 
the data by current processors is often not digestible 
to users. An enzymatic user medol must be developed to 
iaakc the system ccmplete- 

To discuss these points In more detail, first 
refer to Figure 1, which is a familiar block diagram 
of a typical earth resources survey system. The 
system begins with a sensor viewing phenomena of the 
terrain or ocean. If the sensor is a nul tispcctral 
Bconncr, the oh:>crvc.d phens-nena arc the reflected or 
emitted spectral radiance of the scene. Sensor data 
Is collected from reiwtc platforms (hc.nc^ -the nfinui 
remote sensing), so some method of data storage and 
tclcfft^try are required. The next step In the system 
is gGotRc-tric and radiometric preprocessing to permit 
production of a Mp-like reftcfltlmo of tliC true 

of the ievratn. The goal Is to reduce tbc 
distortions gfioftetry and of radiometric fidelity 
to nesllglbl^ VcvcU for further processing and 
analysis. Xhia step in the ptocedwre nty rc<fulTc 
anelliavy (e.gr, knowledge of spacecraft 
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FIGURE 1. TYPICAL EARTH RESOURCES SYSTEM 


attitude and Rround "control points" for performing 
geometric connections and knowledge of atniosphcrlc 
properties for connection of radiometric errors); 
The accuracy and scanner sensitivity to which these 
corrections must be performed are dependent on the 
application being addressed, and may vary widely. 


advancements in extractive processing and preprocessing 
and in user model development which we feci begin to 
close the gap between the sensor’s abilities to collect 
data and the user's desires to digest It, and the 
ability of extractive processing systems to keep up 
with the data. 


The next step in the system Is extractive proces- 
sing, designed to extract information pertinent to 
the user from the data. For example the user may wish 
to know how many lakes there are in Minnesota (claims 
of 10,000 not \v^ichstanding) , and extractive processing 
could be used to map all water bodies and count them. 

The next step in the system is the user model. It 
is a crucial step because it translates the information 
extracted from the remote sensing data, adds ancillary 
informatiorr that a user may have at his disposal, and 
’ creates a product which can help a user make a 
decision. For example, suppose the previously 
tnentloned user wanted to estimate the migratory water- 
fowl produced in Minnesota, Because these birds nest 
in or near water, the amount of water present is an i 

Important but not the only variable. The number of 
water fowl present to breed is important, as are the 
food supply and predator information. These ancillary 
variables arc required to calculate the number of 
migratory waterfowl. . 

Carrying this example a bit further, the user may 
want to know the number of waterfowl to be able to set 
hunting limits In the fall. He will need to consider 
other ancillary information, such as what neighboring 
states arc doing, before making a decision. This is 
the user's decision model. 

do not want to convey the Impression that user 
and decision models are necessarily formal, mathemati- 
cal models implemented on computers. Rather these 
models arc more nearly veil defined procedures that 
managers follow to arrive at conclusions or to convert' 
earth resources processed data to a form they can use. 
The trend is to increasing formalism and mathematics 
in user models v/hcrc inputs and outputs arc usually 
quantifiable. Decision models, with their typical 
socio-economic ancillary inputs, probably will seldom 
be formalized to the stage where computer implementa- 
tion is feasible. 

Vlth this perspective on an oarth resources 
fiystcui, it is perhaps pertinent to point out that 
with increasing developRept of high data rate sensors 
(by HAS/V) and of user interest in outputs from the 
i;ystcri In more tir.’cly fashion (to have maximum impact 
on decisions, inioniation must be timely), the squeeze 
Is on the middle of the system (preprocessing, 
cxtracti\'c processing, and uivor model areas,) In 
the remainder of tlie paper we consider some recent 
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Detection and discriftination of an object by 
multispectral sensing requires differences in the 
radiation received from the object and its surroundings. 
This radiation ’’contrast” is due to differences in 
reflectance, cmittance, or teinperature between the 
object and Its surroundings. 

Dis crimination means devising a decision rule, 
based on measurements from a sample from each of two 
or more given classes, which will enable us to assign 
new samples to the correct classes when we do not know 
to which they belong. Classification means to assign 
samples into groups which sht^ll be as distinct as 
possible. In discrimination the existance of the 
classes Is given; in classification it is a natter to 
be determined. Supervised learning or discrimination 
should achieve greater efficiency because it takes 
advantage of available human knowledge and intelligence. 

In some cases, such human assistance is not possible, 
and an unsupervised learning approach may be required. 

The user defines which process is employed. Having 
made the distinction we use the terms interchangeably. 

The terms identification and recognition will also he 
employed for convenience. 

Basic to this process of discrimination is the 
concept of a signature. In general, a signature is any 
collection of observable features of a material or its 
condition that can be used for precise cla.ssi f ication . 

The features that make up a signature may all be 
obseiTved simultaneously or in a sequence of observa- 
tions spread over a considerable time period. 

A basic element of spectral information extraction 
is the realization that spectral signatures cannot be 
completely deterministic. That Is spectral reflecti- 
vity and cuisslvity measurements of natural objects 
exhibit some dispersion around a mean value (i,c., 
spcctraT signatures arc statistical in character). 

Thus, as we will use the tern, a spectral signature Is 
a probability density function (or set of such functions) 
which characterize the statistical attributes of a 
finite 'set of observations of a material and can be 
used to classify the material or its condition to some 
degree of fineness. 

At tho basis of discrimination theory Is the 
necessity to realize that op Li mum disc riiuinat ion 
techniques require not only that the procedures be 
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tailored to rccop.nlze the Item or material of Interest, 
but also simultaneously, that they be tailored to 
reject other items or natcrials that lie In the vicinity 
of the desired ratcriais but that arc not of Interest, 
the backgrounds In which the items of interest 
arc enbedded. Two types of error are possible: 
failure to claiislfy all of the desired class actually 
present as that class and misclasstfi cation of other 
classes as that class. Photo interpreters commonly 
call these errors of ociission and commission, respec- 
tively. 

Examination of r.easuretnents of the spectral 
reflectivity and omissivity of materials can aief in 
the dc^'clopr.ent of effective discrimination procedures 
by providing insight into the basic optical properties 
of the materials of incerst and their natural varl- 
ability. 

The key to r.ultispcctral recognition is Invariance » 
for example, it is desirable that the classification 
assigned to an object or pattern of interest be inde- 
pendent of the position of that object in the field of 
view, the aspect at which it is viewed; the background 
against which it is seen, partial obscuration of the 
object, minor changes within a class, and changes in 
illumination or atmospheric condition. It is not too 
difficult to provide any one of these invariances. To 
provide all of the desired invariatices with a practi- 
cal amount of hardware, however, requires that the 
preprocessing and feature extraction mechanism extract 
the essence of the classes to be identified. 

In a multlspcctral recognition system, five major 
functional divisions must be considered. The input 
picture element or pixel is a vector quantity made up 
of many components hence the dimensionality of the 
Input space may be large. 

The purpose of signal conditioning, or **preproces- 
sing” Is' to provide data preparation and handling, 
to provide geometric and radiometric correction, to 
provide a convenient input format, to provide invari- • 
ance , to provide in nany cases a reduction in the 
ditnensionali ty of the input data, and most importantly, 
to emphasize or enhance aspects of the input signal 
vhlch are deemed important. 

An almost universal approach to recognition Is to 
extract properties or features from the original 
signal, and to perform the recognition on the feature 
profile of the input signal. This serves several 
functions. First, by reducing the input pattern to its 
essential features, the memory required for storing 
the signatures is reduced. Secondly, by reducing the 
input pattern to independent features, a considerable 
Binount of invariance to exact torn Is obtained. 

Finally* a degree of invariance to noise and background 
tiay be achieved. 

Host decision mechanisms are based on mul tivariate 
discriminate analysis that partitions measurement 
cpacc on the basis of the training set signatures which 
then allow's a decision to be made for an appropriate 
classification- for each input pixel. 

The function of dij;play is attendant with each and 
every one of the other functions of a multispcccrnl 
processing sysicm as one facet of 'the man-machluc 
Interaction. Interactive controls and co»n»ands allow 
the intervention of the operator to direct that certain 
t!^^lXgs be done which w'ould not otherwise be done auto- 
tntically. 


The Multivariate Interactive Digital Analysis System 
(MIDAS) system is an attempt to solve the pr*.'blcrn 
of real’ time multlspcctral data processing in an opera- 
tional system. We have been aware of this problem since 
the initial design and operation of the H-7 scaancr and 
processing equipment at ERIM over ten years ago, A 
rcal-tinc processor (S?^\RC - SPcctral Analysis and 
Recognition Computer) was completed in 1967, allowing 
parallel decision operations using a nultlvariatc maxi- 
mum likelihood algorithm thus making possible for the 
first time, rapid processing of 12-band data (Ref, 1}. 

It was evident that a considerable amount of assistance 
and cine was needed by the user in setting up such a 
machine and controlling its operation. This led to 
the evaluation of a hybrid system employing a general- 
purpose digital computer for control by the user of the 
aystem (Ref. 2). Processing of remotely sensed data is 
an interactive process in which the man and machine 
Twst , in fact, be considered as the real processing 
system. 

This is not apparent In general purpose computer 
processing systems using software classifiers since the 
machine is so slow by comparison to parallel 
approaches that an operator ic easily able to keep 
with Che task (hours or days is the normal time intervnl^ . 
In this MIDAS system, where the parallel-pip«?Unc 
processor is substantially faster, the time I'cqx^ired by 
an operator is three or four times longer than chat 
used in processing. Increases in processor speed wi li- 
chen provide little improvement in chrougbpu*:. Tr- 
becomes evident that yell designed, interactive display 
and control subsystems will, in reality, offer the 
greatest gains in throughput. 

Froni another point of view, the objective of 
providing greater speed also provides much lower cost. 
This is true directly for the operational situaticn 
but also for research and development. It has been 
estimated that processing costs in an operational 
system can. be reduced by about a factor of tv;enty or 
more from some present day feasibility processing 
costs. ^ ♦ 

The overall system hardware is sho™ in Figure 2. 

The MIDAS system (ref. 3) consists of several princi- 
pal subsystems; the general purpose computer (DEC-PDF- 
il/45), used for analysis and control but not preproces- 
sing or classification, the classifier, the prcproces*^ 
sor, the Input subsystem, the control subsystem, and 
the output displays. Of these, the PD?-ll-/^5 computer, 
the classifier, the control system, and part of the 
input system have been designed and tested in the 
first year of the NASA A.XFE supported program. The 
remaining units to. be added are indicated as cross 
hatched blocks in Figure 2, 

The Input subsystem' inputs data stored on 1) High 
Density Digital Tape (BDDT) , 2) Computer Ccnipatible 
Tape CCCX) , and 3) Analog Tape. The HDDT input is the 
principal input modi uni because of its high data rate 
and digital format. Data stored on HDDT^s at 20,000 
bits per inch anciplayed back at 120 inches per second 
gives a data rate of 2.4 megabits per second channel. 

For ten channel data this is a bit rate of 24 ridigabits 
per second. Another advantage of UDDT is the efficient 
storage of data and a KASA st-ondard format for Earth 
Resources Surveys has been ‘ rccpntf?5Cndod. The data 
contained on an HDD! nay be as much as on fifty CCT's. 

The second input medium is a 7 or 9 track multi- 
density CCT, The. availability data from many 
sensors on a common format is desirable, however, the 
tJnc to produce such tapes from tl;c original data b63 
been quite long and may tiioke the data out of date. 
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The data rate froin CCT Is at least an order of magni- 
tude below the systcn capacity. A 58 megabyte disk is 
provided on the PDP-ll/4i5 for data storage and storage 
of classification results. 

The third medium of data input is from analog tape 
through the use of multiple A/D converters. Although 
this is an efficient data. storage medium capable of 
xnaxlmum data race for the NIDAS system, it appears as 
the third choice bec.^use few MSS systems record data 
in this form. 

The preprocessor performs signature extension 
preprocessing calculations as needed. These are; 

1) calculation of scan angle correction functions, 
both additive and multiplicative; 2) selecting and 
calculating channel ratio transforins, and 3) calcula- 
ting a linear dimension reduction transform. 


bility for the classes of interest to the user. This 
has the desirable effect that the classifier can per- 
form a classification operation in which the accuracy 
:of classification using the smaller number of dimen- 
sions is equivalent to that obtained with a larger 

number of untransformed dimensions. 

> 

The classifier performs a maximum:- likelihood 
decision, assuming a multi-model Gaussian multivariate 
distribution. This assumption has been well justified 
at this time by over 100 experiments using multispcctral 
data at ERIM (ref, 2) and, as time goes on, by more and 
more experience at NASA and other centers. Although 
‘simpler algorithms can perform well for some data 
sets, a significant percentage of applications demand 
this po\t?erful decision rule. No penalty in speed and 
only a small additional cost occurs for using this 
algorithm. 


The calculation of scan-angle correction functions 
is a step normally required (ref. A) to remove path 
radiance or other atmospheric and angle dependent 
variations. 

Ratio t rans forms (ref. 5) may be used in pre- 
processing to provide data which is relatively Insensi- 
tive to illumination and transmission variations, where 
transference of signatures from one frame to another 
Is desired, or v:hcrc spectral features can be enhanced. 
There is no clear-cut means of examining a data set 
to decide, in’ a priori manner, v)iich of several possible 
transforms is needed for a particular scene. The 
inctbod used would norr,:ally be one in tvhich each of 
several transforms arc performed oi\ the training sets 
and tlie resulting data is tested for the optimal, 
probability of correct classif 1 cation using the 
training sec and test set data to choose the transform 
to* be employed. 

The pui'fiOSC of tlic linear transform (ref. 0) Is 
to provide a new set of data in which the spectral dat«a 
is conbincd in sucli a manner that the d i.mcns ion- reduced , 
transformed data has essentially tlic same discrlmlna- 

/♦ 


The actual computations are performed by sets of 
time-shared arithmetic units arranged in parallel 
pipelines. Each stage supplies its results to a 
subsequent stage for further processing. Precision 
varies between 8 and 16 bits ns the data progresses 
through the pipeline and acquires greater significance. 

• It^will accept 8 input signals and can classify 
the results into one of 9 classes (including the null 
class) at its output on one node of operation. In 
another mddc of operation it can accept four ifiput 
signals and classify the results into one of 17 ci isses 
at its output. The decision rate Is 200 x lOVsecond 
allowing classif ictation of an entire ERTS frane in 
less than one minute. 

A three color video display using hiph resolution 
a shadow-mask CRT and n MOS storage for display- 
refresh offers an interactive capability. The unit 
chosen .at this tine Is the GX- 100/ 2 00 series 

configured to hold n 5 bit color vector with 512 x 512 
elements on the display. 
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The displAy offers Iraacc display, alphanutncrlcs , 
vector generation, 2 inovablc cursors, digital zoom 
and wove or moving window display and table loo ~“P 
of predetermined colors. This color display prov es 
the major Interface for a fast man-machine interface. 

A eccnc or portions of a scene can be displayc 
alEupst iirjtediatciy from disc as the user requests t* 
By havtufi the ability co chance colors quickly on the 
CRT throuB^^ of the table look-up » enhanced 

pictures can be viewed which help locate training 
sets of arbitrary shape. These, enhanced pictures 
laay also be put in hard copy form* 

A second interactive CRT Is part of the system for 
alphanumeric information display. 


Socctral-Snatlal Processing 

ERIM has explored both applications of Inforoatlon 
from ERTS, with pronislng results. In one o 
V. h.'v. .r..d 

the kev Step Is the formation of spatial 

l!; .ttrlbuce. .1 .h. FoU~.bE .he 

Kansas, spatial features T^cse 

over Michigan .ns shoi/n in Table 1. (Ref. /)• 

spatial features were formed by measuring 

ERTS band MSS-7 signal level in a 9 x 

the pixel of interest at the center. Then 

specLol and spatial siganturcs „ 

tLraln categories in the Ann Arbor-Brighton, Micnigaa 


An extremely important device for the MIDAS system 
is a fast color printer. There have been a number of 
methods proposed for this function Including CRT-color 
filter camera devices, multi-laser film devices, LED 
color array systems for printing on film and colored- 
'Ink printers. Of these, the most attractive method 
Is the ink-jet printer since it produces a usable out- 
put immediately with no delay for color processing. 
This appears capable of producing a picture of about 
10^ elements in 10 minutes or less. An alternative 
possibility is a multi-LED scanning film printer 
using a Polaroid or Eastman fast color development 
method. ' ' 


Having this fast turnaround for intermediate or 
final recognition results will pennit the user to be 
much more efficient. This is especially true when 
enhanced false color maps of training, areas can be 
made quickly. Present ERTS data sets can take up to 
40 tnan hours or analysis to locate training sets 
because of low contrast on single channel black and 
vhite hardcopy, and because of the poor quality of 
available hardcopy. 


We next ran an optlinum feature selection proptcim 
based on an .algorithm which selects ’ 

along with the channels .already seiecccu, . 

the average pain-.-isc probability of miscias^ificaui 
between pairs of sigt^atufes. 


The results of the 


optlf-um channel ordering are presented in T^blc 2. 
Sotc that 2 of die first 4 fe.aturcs arc spatial .e. 


eaturcB. 


TABLE 1. SVATIAL FEATURES DERIVED 

1-4 mean signal level in each ERTS channel over 
9x9 array {P4 - wy) 


'5-8 standard deviation over 9x9 array in each 
ERTS channel (04 - 07) 

9-14 normalized indicator of Intet-channel covariance 
between ERTS channels over 9x9 array CR45“^67) 

15-22 average power in 4 spatlk frequency bands (from 
1/2 pixel to 1/5 pixel) for bands MSS-5 and 
MSS -7 CQ25- Q55. snd “ ^57^ 


In developing a complex system of digital hardware 
- such as MIDAS, it is necessary to provide a method or 
Eet of methods which provide a quantitative report on, 
the performance of the hardware for fast, low cost 
maintenance. This set of methods comes under the 
general category of machine diagnostics. 

A substantial number of diagnostic programs have 
been developed which make use of the diagnostic bus 
designed into the system. The set of these diagnostics 
can be broken dowm into the following four types: 

1) an automatic in-piace fault isolation diagnostic, 

2) classifier sub- assembly diagnostics, 3) off-line 
individual card diagnostics, and 4) a classification 
accuracy analysis package. 

recent PE Vi:Lui'MEN'T5i IN PROCESSING METUODOIOGV 


table 2. ORDERING OF SPATIAL-SPECTRAL FEATURES FOR 
terrain classification using JUNE MICHIG 2LN 


ERTS DATA 


of Channels 

Feature No* 

Prob. of Misclass 

1 

1 

77 

2 

11 

41 

3 

20 

24 

4 

4 

18 

5 

8 

14 

6 

2 

13 

7 

22 

12 

B 

• 10 

11 

9 

5 

10 

10 

6 

8 

11 

3 

8 

12 

9 

S 


The launch and subsequent long operation of the 
ERTS-I nul tlspecrfaJ scanner have allowed users to 
ohs-^rve areas of interest rcpentedly at 80 meter 
resolution. Because the ERTS MSS has four rather 
broad spectral bands, extractive processing using 
only spectral ion fi^quently yields ifnperfcct 

GCpnration of classes o£ materials of interest. 
Partially to oVercoRd teKc performance ERTS 

crcctraL channeie relative to those available form 
typical aircraft scanners, users have exploited the 
spatial information inherent in the BRTS data along 
witi^ the spectral data. Further, using the repetitive- 
coverage capability of F.RTS, some users have used the 
temporal' variation of cpcc-tral data as inputs to the 
pattern rccosnitlcn processors. 


*See Table 1 for identify of features. 

Table 3 sumrvirlzcs the pcr£or«\ance of the pattern 
recognition classifier using tht twenty-two features 
of Table 1. From the nuinbcr of spatial features used 
in the classification, their position in the list of 
optimum channels, and the performance of the classifier 
using the spatial features, one nay obtain some 
intu/tion of the utility of spatial feature?, when used 
with spectral features in pattern recognition, 

VTillc it would be premature to suggest that this 
type of fipatiad-'Spectra^ processing we hava perfotned 
using ER.TS data t-ilU be cltlrtately useful in oU 
discipline dtedd particularly dac to its coarsening or 
the resolution, we still cone loiia that the combined 
use of EpatlaU and spectral feeturce In conventional 
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TADLE 3,. I’ERFORi'LANCE OK CUSSIl’IER USING SPECTRAL/ 
SPATIAL KEATURES 
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pattern rccocnition algarithias has definite promise “ . 
and should be explored further. The work described 
here was pcrfor;T^cd for the U. S, Army-MERDC, under . 
subcontract to ERIM from batellc Memorial Labs* 

Spectral»-Tor.poral Processing ' 

Because of the periodic coverage of ERTS, the 
temporal variation of spectral signatures can be 
utilized to help discriminate objects, ^ This capability 
Is enhanced with ERTS data over the previously available 
aircraft data because the periodic coverage is obtained 
from a sun**synchronous orbit. Further » the data arc 
obtained from a relatively stable platform and at small 
scan angles from the nadir (typically + 5,5°). Both 
facts greatly ease the data registration problem •— 
that of merging data collected at two different tiines 
so that pixels arc precisely overlaid. 

To illustrate the advantage of rnultitcmporal 
spectral processing, consider the following example 
of mapping of natural vegetation in northern Michigan. 
(The work being discussed here was performed for KASA 
under grant NGR23-005-552) , 

Tabic A shows the signatures of a mixture of 
hardwoods, conifers, and grass and shrub swamp on two 
different dates — June and March, Shot-m are the 
mean values of the signatures, with the standard 
deviations in parenthesis. Notice that in June, the 
two classes signatures overlap appreciably in each 
ERTS channel -- the mean difference between 
classes is less than the standard deviation in each 
ctiannel, A pattern recognition device hopelessly 
confuses these tv^o classes in June. 

But in Marcli, the two classes are more separable, 
as shown in tlie second half of Tabic A. This 
separation of scjtic classes at one time of year and 
not at others can be exploited to improve maps of 
natural vagfetatioa areas. the challenge is to un- 
derstand Vegetation &igt^atvres diangc wiLl» time, 

BO that VC may specify exactly when data arc to be 
collected. 

A second advantage of the periodic coverage Is 
the ability to monitor the development of agricultural 
crops. While coverage over agricultural 

areas significantly influences the amount of periodic 
coverage available, the cov&rdge can still be used, 

T)ie potential exists for crop yield osti [nation and 
prediction, c^f^p&rlng the crop development, as sensed 
from EKIS, with a crop calendar normal devc Jopmetf^t. 

can accept either spectraH-fipatlftl or spcctrnl- 
tcmporftl data BS easily as sped re I alone. 


TABLE 4, SIGNATURES OF TWO VEGETATION CLASSES ON TWO 
DAYS. 


June 

MSS-5 

HSS-5 

MSS-6 

HSS-7 

l!dw/ Coni f /eras s 

27.32 

(0.91) 

17.82 

(1.45) 

49.62 

(5.52) 

27.53 

(3.6.8) 

Shrub /Swamp 
March 

27.69 

(1.20) 

17,25 

(1.53) 

47.06 

(5.69) 

25.75 

(3.94) 

Hdw/Conl f /Grass 

29.96 

(3.75) 

28,07 

(3,11) 

29.93 

(2.96) 

15.82 

(1.28) 

Shrub/Swamp 

23.69 

(1.74) 

20.25 

(4.04) 

22.19 

(3.60) 

11.38 
(2.60) . 


Mean signature values are shown, with standard 
deviations in parenthesis. 


USER MODEL DEVELOPlfENT 

One critical portion of the Earth Resources 
Survey System is the User Model, which relates Che 
output of the extractive processors and ancillary data 
to generate information which a user can employ directly. 
User models may be very simple ~ if the user wants a 
laap of vegetation types, the output of the extractive 
processing may directly serve him directly, and the 
user model is absent. But If a Department of Agri- 
culture official wants to know what is the projected 
wheat production in Kansas, Che user model may 
combine the total productive acreage of wheat (obtained) 
from a remote sensing system), with some farmers^s 
estimates of the yield of particular fields, and some 
estimates from the weather service of future weather 
trends, to calculate the total production of wheat in 
Kansas. 


As an example of a user model, consider 
Figure 3. This model predicts the population of 
migratory waterfowl, given the water supply conditions, 
the food supply conditions, and ancillary variables 
such as number of nesting pairs, predation, and 
mortality. Also shexim in Figure 3 on the right are 
elements of the decision model used by the USDI Bureau 
of Sports Fisheries and Wildlife, The model shown 
is, at present, conceptual and represents the work of 
BSFl^ personnel, ERIM and BSR^ are working jointly 
(under NASA funding from ERTS and Skylab programs) to 
develop remote sensing inputs to this model. 

Figure 3 shows at the left, inputs to the user 
model. Asterisks indicate those types of information 
obtained from low altitude aircraf t-bomc observers 
(supplemented by grouTid observations), and data 
potentially derivable from remote sensing techniques. 
The goal of the user model Is to estimate the fall 
population of wature *and i^m&turc birds. Old birds 
arc estimated by low altitude aerial survey and 
ground counts of nesting pairs on a sampled basis. 
Estimates of su«mcr r*o reality are also inade from 
ground observations and past experience. The number 
of breeding pairs, coupled with the May ond July pond 
numbers arc to estimate the number of new docks. 

This info relation is augifjentcd by the number of broods 
obtained on a campled basis from ground survey, 

HuntJng rer.ul.it ions arc defined based on the 
pop ul.it Ion size, the of J^arvestlng ducks 

in Can^dd, ond the tarrying capociLy of the habitat- 
Remote sensing has cn Impact In assessing capacity of 


G 
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habitat, especially In the assessnont of quantities 
end distribution of natural vegetation* 

As nn exontplc of the annual production equation. 

Figure 4 shows n linear equation in several variables 
vbich BSFW personnel have derived to predict new 
production of r.allard ducks* The coefficients of 
the raodel are derived on the basis of historical 
experience* If similar models can be constructed 
for other species, then Che general user-decision 
model of Figure 4 can become more operational. 
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DATA INPUTS 
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FIGURE 3 

DETERMINING HUNTING REGULATIONS BASED UPON THE ESTIMATED 

magnitude of the fall flight of migratory waterfowl 

♦Obtained from aerial observations and ' • 

adjusted based on scltclod groui'.d counts 

♦♦Obtained from aerial observations -^data ‘ ■ 

potentially derived using remote sensing ’ , ' 

techniques 


Y = 7.920 + 1.468 Xj - 0.624X2 - O.O 28 X 3 + 0.016 X^ 

where Y = Predicted number of mallard young (millions) 

Xj = July ponds (millions) 

Xg = continental mallard breeding population (millions) 
Xg e percent of ponds remaining from May to July 
X.» - index io number of broods (Hi ou sand 5; unadjustci^ 

EXAMPLE OF A MODEL FOR PREDICTING 
ANNUAL PRODUCTION OF YOUNG MALLARDS* 

♦After Ccis, A. D., R. K. Martinson, and D, R. Anderson (1969) 


FIGURE 4 
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